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Abstract 

We summarize our current understanding of the connection between the QCD phase 
line and the chemical freeze-out curve as deduced from thermal analyses of yields of 
particles produced in central collisions between relativistic nuclei. 

Quantum chromodynamics (QCD), the theory of strong interactions, contains, as a key 
prediction, a new phase of matter, the quark-gluon plasma (QGP). This partonic QGP 
was the main ingredient (in addition to leptons and neutrinos) of the matter in the early 
universe until about 10 /is after the big bang. Around that time the QGP was transformed 
into hadronic matter through the QCD phase transition. In general, strongly interacting 
matter in equilibrium is characterized by two quantities, the temperature T and the baryon 
chemical potential /i^ or equivalently, the (net) baryon density Ub- 

How does this phase transition come about? Due to asymptotic freedom, the strong 
coupling constant a s runs, i.e. diminishes with increasing energy scale, implying that in- 
teractions among strongly interacting particles will get weaker as T or \xb (or both) are 
increased. One may ask at what energy scale will there be new physics such as deconfine- 
ment for strongly interacting matter in equilibrium? The only two scales relevant appear 
to be Aqcd ~ 200 MeV along the temperature axis and the nucleon mass along the 
fiB axis. From phenomenological considerations and based on ever more accurate solutions 
of QCD on a discrete space-time lattice ('lattice-QCD') convincing arguments were put for- 
ward to demonstrate that strongly interacting matter undergoes indeed a (phase) transition 
from a dense hadronic medium where all constituents (hadrons) are confined to a decon- 
fined plasma of interacting quarks and gluons. This implies the existence of a line in the 
(T — ub) plane, the QCD phase boundary, anchored by parametrically critical parameters 
T c ({ib = 0) ~ Aqcd and (//b) c (T = 0) ~ m^. For recent detailed reviews on the QCD phase 
diagram, where many of the above short arguments are exposed in detail and the state of 
theory and relevant experiments is summarized, see [TJ [2] . 



This article is an attempt to update a paper on the subject which we wrote 15 years ago 
[3] in honor of Gerry Brown's 70th birthday. While we have (with one exception [I]) never 
published together with Gerry, his constant probing and provocative questions on the QGP 
and related areas have, during nearly 2 decades of close scientific interactions, deeply influ- 
enced our thinking on the subject. The experimental basis for making a connection between 
data from ultra-relativistic nuclear collisions and the QCD phase boundary is dramatically 
improved after 15 years of intense research on the subject. Given all this we consider it an 
honor and an opportunity to be able to summarize our most recent thinking on the subject 
as part of a Festschrift on the occasion of Gerry's 85th birthday. 

Strongly interacting matter at high temperature and baryon density is produced in col- 
lisions between atomic nuclei at ultra-relativistic energies. Depending on the center-of-mass 
energy reached in the collisions the temperature and baryon density can be 'tuned' system- 
atically, as discussed in detail below. A sequence of such experiments conducted over the 
past 25 years over a wide range of energies has provided convincing evidence for a new state 
of matter, see, e.g., [3 El El El El HD1 HH H2]- In November 2010, first results from Pb-Pb 
collisions at the Large Hadron Collider at CERN have provided a first glimpse into fireballs 
with ultra-high temperatures (approaching 1 GeV) and opened a new era in quark matter 
research [Ml M> HHJ UM El EE] . 

The range of temperatures and densities reached in these experiments is clearly of the 
order of or (for RHIC and LHC experiments) significantly exceeding the typical values dis- 
cussed above. For a review and summary of the evidence, see [19]. Much harder has been 
to assess the degree of equilibration reached in the collisions and to provide direct evidence 
on the QCD phase boundary. Below we will argue that the most direct information on 
equilibration and the QCD phase boundary available to-date is obtained from analyses of 
the multiplicities of hadrons produced in central collisions between ultra-relativistic nuclei. 

These data have been analyzed independently by various groups, using hadron resonance 
gas models with chemical freeze-out as a key ingredient. In this resonance gas approach, 
described in detail in [20], it is assumed that the final hadron yields result from the decay 
of fully equilibrated hadronic matter, comprising the full QCD hadron mass spectrum - and 
thereby implicitely containing all of QCD. The measured hadrons, those not decaying further 
by the strong interaction, are then obtained by summing up all strong decay products of 
the hadron mass spectrum with their corresponding thermal weight. By assuming that the 
measured hadrons are generated at a common surface at which all particles simultaneously 
decouple, values of the baryon chemical potential, and temperature, T, on this surface, 
the chemical freeze-out surface, are extracted. Fitting these two parameters, fib and T, 
together with the volume parameter gives values for the particle abundances which are in 
very good agreement with experiment [2H [221 E31 EH [251 EHl EH EU [291 EQl EH E21 [331 EU, 
for an overall summary see [2D] . 

These results clearly demonstrate that, at chemical freeze-out, a very high degree of 
equilibration is reached in central nucleus-nucleus collisions, implying that matter in the 
proper sense with thermodynamic properties has been formed at the latest at this stage of 
the collision. We note that the system most likely is in or very close to equilibrium much 
earlier in the collision, as demonstrated by, among other things, hydrodynamic analysis of 
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Figure 1: The temperature and baryon chemical potential of Statistical Model fits to hadro- 
chemical abundances as a function of center of mass energy per nucleon pair for collisions of 
heavy nuclei (Figure taken from [33| 134"]). 

azimuthal anisotropics ('elliptic flow'). For a brief summary see |12[ [19]. 

The resulting values of \ib and T are shown in Fig. [1] as functions of center-of-mass energy 
per nucleon pair. We note first a very important outcome of these investigations: near 10 GeV 
center of mass energy, the temperature saturates with increasing beam energy, reaching an 
asymptotic value of about 160 MeV, while the baryon chemical potential decreases smoothly. 
This temperature saturation concept can be tested thoroughly when the newly taken data 
at the much higher energy of a/stvtv = 2.76 TeV from the LHC are analyzed. 

All these data points have been obtained by analyzing hadron yields, implying the exis- 
tence of a limiting temperature to which a hadron 'gas' can be heated, as has been conjec- 
tured more than 40 years ago by Hagedorn [35]. The constancy of T above ^snn = 10 GeV 
indicates that a boundary has been reached. Is this the QCD phase boundary? Could there 
be hadronic matter hotter than is observed here? Since the data analyzed here comprise the 
full hadronic spectrum, this would imply the existence of a hot and dense hadronic medium 
between the QCD phase boundary and the chemical freeze-out curve that leaves no trace 
in hadronic observables. As we will discuss below, this would contradict the existence of 
chemical freeze-out outlined above. 

Plotting these temperature-chemical potential pairs for all available energies results in a 
phase diagram-like picture as is illustrated in Fig. [2J In the \xb region from 800 to 400 MeV, 
as T increases from 50 to 150 MeV, the experimental points rise approximately linearly. In 
contrast, below fis — 400 MeV, the temperature is approximately constant, T ~ 160 MeV. 
The highest collision energies studied to date at RHIC are those for which \ib — 25 MeV. 
Soon an experimental point at \Ib — 1 MeV will be available from the recent campaign at 
the LHC. Also shown on this plot are lines of fixed energy per particle and fixed entropy 
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Figure 2: The decoupling temperatures and chemical potentials extracted by Statistical 
Model fits to experimental data. The freeze-out points are from Refs. [28] and (33J [36j [37] . 
The open points are obtained from fits to data at mid-rapidity, the full-points refer to An 
data. The inverse triangle at T = indicates the position of normal nuclear matter at T = 0. 
The lines are different model calculations to provide a phenomenological understanding of 
the freeze-out curve [32J ESI EH] • The shaded lines are drawn to indicate different regimes in 
this diagram. Figure taken from [ID] 

density per T 3 as well as a line of hadron percolation (32j [38j [39] . 

These experimental results can be compared to the phase boundary computed on the 
lattice [HI H2]. Numerical simulations in lattice QCD can be performed at nonzero tempe- 
rature only for small values of fis, typically fis < T, without running into the 'sign' problem. 
At fiB = 0, these simulations indicate that there is no true phase transition from hadronic 
matter to a quark-gluon plasma, but rather a very rapid rise in the energy density at a 
temperature T c of 155 — 175 MeV within current systematic errors, in very close agreement 
with the chemical freeze-out temperatures determined for small values of \ib- 

Further, studies using the lattice technique imply that T c decreases very little as hb 
increases, at least for moderate values of hb- A very recent study [43J provides a quantitative 
estimate of the curvature of the phase boundary near \i B = 0. By comparison to the curvature 
of the line of fixed energy/particle [3S] which has been constructed to get a phenomenological 
understanding of the complete chemical freeze-out curve, these authors argue that there is 
a significant difference between the data on chemical freeze-out and the lattice predictions, 
implying in their view the existence of a dense hadronic phase between chemical freeze-out 
and the QCD phase boundary. 
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A precise determination of the curvature of the chemical freeze-out curve can, however, 
only be obtained by analyzing the data themselves, not a curve such as the line of fixed 
energy/particle or any of the other lines shown in Fig. El We have recently investigated 
jH] the curvature of the chemical freeze-out curve in the region \xb < 250 MeV. The results 
are presented in Fig. [3J The chemical freeze-out points are taken from [ID] . In addition, we 
show in Fig. [3J also fit results from the STAR collaboration [25], where only pions, kaons 
and protons are analyzed to determine the freeze-out curve. This leads to a slightly lower 
chemical freeze-out temperature, but we believe that its energy (or fis) dependence contains 
the correct information. Independent of which chemical freeze-out parameters are used, the 
resulting description using the parabolic fit function as in [4_3] (long dashes) and the same 
function with twice the slope parameter (short dashes) are in agreement with the data, 
within the (still significant systematic) uncertainties, as shown in Fig. [3j For reference we 
include, in Fig.|3j the line of constant E/N [39J. From the available data one cannot construct 
evidence for a dense hadronic phase between chemical freeze-out and the phase boundary. 
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Figure 3: Comparison of the experimental /Ib dependence of the chemical freeze-out tem- 
perature with recent lattice predictions [13] . The STAR chemical freeze-out values are from 
[15], the other values are from [lUJ. For more discussion see text. 

With the parametrizations of T and \i& from Fig. [T] one can compute the energy de- 
pendence of the production yields of various hadrons relative to pions, shown in Fig. HJ 
Important for our purposes is the observation that there are peaks in the abundances of 
strange to non-strange particles at center of mass energies near 10 GeV, i.e. where the tem- 
perature reaches its limiting value. In particular, the K + /tt + and A/ir ratios exhibit rather 
pronounced maxima there. We further note that in the region near 10 GeV, there is also 
a minimum in the chemical freeze-out volume [3TJ H6] obtained from the Statistical Model 
fit to particle yields [3TJ [33] , as well as in the thermal freeze-out volume obtained from the 
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Hanbury-Brown and Twiss (HBT) radii of the fireball [37]. The energy dependence of the 
chemical freeze-out volume is shown in Fig. |5j Included in this figure is the most recent 
point at LHC energy from ALICE [13] obtained by analyzing the pseudo-rapidity density of 
charged particles for central Pb-Pb collisions. 
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Figure 4: Energy dependence of hadron yields relative to pions. The points are experimental 
data from various experiments. Lines are results of the Statistical Model calculations. The 
Figure is taken from (33J I34"]). 

These experimental observations have long resisted interpretation in terms of a transition 
between hadronic matter and a quark-gluon plasmsQ- The general structures observed in the 
data are well reproduced only by the most recent model calculations [33] • There, it is shown 
that these structures arise due to the interplay between the limit in hadronic temperature 
(see Fig. [1]), in our interpretation due to the QCD phase transition, and the rapid decrease 
of fiB with increasing energy. The minimum near a center of mass energy of 10 GeV and the 

■"■We note the interpretation given in [48] , obtained within a schematic 1st order phase transition model. 
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rise towards higher energies in the thermal freeze-out volume jU] was explained as due to the 
increasing meson to baryon ratio in this region combined with a meson-meson cross section 
which is significantly smaller than the pion-nucleon cross section. The chemical freeze-out 
volume depicted in Fig. [5] behaves similarily, with a minimum at an energy close to where 
the chemical freeze-out temperature saturates. At the highest (LHC) energy the chemical 
freeze-out volume reaches, for Pb-Pb central collisions, a value of about 4 times, the thermal 
freeze-out volume about 7 times that of a Pb nucleus. The strong volume increase between 
RHIC and LHC energy, where neither T nor \xb are expected to change much, provides 
a strong argument that chemical freeze-out takes place at a fixed density. These findings 
further strengthen the connection between hadron gas and quark-gluon plasma. 

The same chemical freeze-out mechanism also governs the production of light nuclei 
and their anti-particles and even complex and exotic objects such as light hypernuclei and 
their anti-particles [4*9| l50| l5Tj . This comes at first glance as a surprize to many, as the 
binding energy of these objects is typically 1-2 orders of magnitude smaller than the chemical 
freeze-out temperature of T « 160 MeV. Can it be that such loosely bound objects are 
messengers of the phase transition between hadronic matter and the quark-gluon plasma? 
Entropy conservation comes to the rescue! After the phase transition the hot fireball expands 
adiabatically and the overall entropy as well as the entropy /baryon are conserved quantities. 
As was realized already more than 30 years ago [52] the assembly yield of complex nuclei from 
a gas of hot nucleons is a measure of the entropy/baryon, and hence can be used to diagnose 
the hot and dense phase of the collision. The answer to the above question is therefore 
strongly affirmative: the yield of composite objects such as light nuclei and hyper-nuclei 
(and their anti-particles) produced in central nucleus-nucleus collisions at ultra-rlativistic 
energy is a direct measure of the entropy/baryon at chemical freeze-out, and hence, of the 
QCD phase transition. 

At RHIC energies, chemical freeze-out was shown [53] to take place very close (within less 
than about 10 MeV) to the phase boundary, driven by multi-particle collisions in the high 
density regime of, and by the the rapid density change across the phase transition. Further it 
is argued that freeze-out ends when the system is fully hadronized, i.e. at comparatively low 
density in the hadronic phase. Without a phase transition, freeze-out in a purely hadronic 
medium takes place over a considerable time and temperature range, as has been recently 
demonstrated from very general considerations [51]. This would necessarily lead to different 
freeze-out parameters for each hadron species due to widely different hadronic cross sections. 
This is not observed. Conversely, the rapid change of density with temperature near the 
phase transition makes the system insensitive to the different cross sections. We conclude 
that the observed (nearly) simultaneous chemical freeze-out of all hadrons is inconsistent 
with hypothesis of a dense hadronic phase between the chemical freeze-out line and the 
phase boundary, at least for relatively small values of \x& < 250 MeV. 

We believe, however, that the above argument is generic [53]: to ensure simultaneous 
freeze-out (within a very small interval in temperature and chemical potential) of all hadrons, 
the freeze-out curve has to be very close to a line with a rapid density change, significantly 
faster than can be obtained within the scenario of an expanding hadronic medium. An 
immediate consequence of this would be that the chemical freeze-out curve delineates phase 



7 



E 

^5000 



I I I I 1 1 1 1 — III 



4000 



3000 



2000 



1000 



central collisions 



i i 



10 



1<f 



10 



Figure 5: Energy dependence of the volume for central nucleus-nucleus collisions. The 
chemical freeze-out volume dV/dy for one unit of rapidity is taken from Ref. [21] • The LHC 
point is from analysis of the recent ALICE data |13j . 



boundaries, not only for small values of \ib but everywhere. 

But what provides the phase boundary for large values of /ie, where the deconfmement 
transition seems far away, at least if one follows the guidance from lattice QCD calculations? 
In line with Gerry Brown's daring attitude concerning new scientific directions we have 
recently speculated [40J that the transition from hadronic to quarkyonic matter might provide 
the missing link. Another, related possibility might be that, at high baryon density, the lines 
corresponding to the chiral and deconfinement phase transition split, leading to a region of 
chirally symmetric but still confined matter in the QCD phase diagram. 

Whatever the case may be, it is clear that the phase diagram at high baryon density 
is difficult to explore, both experimentally and theoretically. Speculations have been put 
forward about the possible existence of a critical endpoint [55] or a triple point [ID] in 
the QCD phase diagram. However, despite significant efforts, no clear signals have been 
discovered to-date [561 EZ] °f either critical points or, more generally, of dense phases beyond 
chemical freeze-out. This research on the high baryon density region of the phase diagram 
remains a challenge for the future and will have to be largely driven by the running campaign 
at RHIC [57J and the planned new experiments at FAIR [58] and DUBNA [5J] . 

We have provided strong evidence that, for not too large values of the baryo-chemical 
potential \xb < 250 MeV, the experimentally observed chemical freeze-out curve closely coin- 
cides with the QCD phase boundary between the hadronic world and the quark-gluon plasma. 
This implies a direct connection between results from experiments with ultra-relativistic nu- 
clear collisions and a fundamental prediction of QCD concerning the phase structure of 
strongly interacting matter. For large values of hb, i-e. large net baryon densities, we can 
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presently only speculate, but in Gerry Brown fashion have put forward an argument that 
also there the chemical freeze-out curve is driven by a phase transition. 
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